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Abstract: The multi-parameter character of supersymmetric dark-matter models implies the 
combination of their experimental studies with astrophysical and cosmological probes. The 
physics of the early Universe provides nontrivial effects of non-equilibrium particles and 
primordial cosmological structures. Primordial black holes (PBHs) are a profound signature 
of such stmctures that may arise as a cosmological consequence of supersymmetric (SUSY) 
models. SUSY-based mechanisms of baryosynthesis can lead to the possibility of antimatter 
domains in a baryon asymmetric Universe. In the context of cosmoparticle physics, which 
studies the fundamental relationship of the micro- and macro-worlds, the development of 
SUSY illustrates the main principles of this approach as the physical basis of the modern 
cosmology provides cross disciplinary tests in physical and astronomical studies. 

Keywords: Cosmology; Particle physics; Cosmoparticle physics; Supersymmetry; 
Primordial Black holes; Non-equilibrium particles; Antimatter 


1. Introduction 

The multi-parameter space of Supersymmetric (SUSY) models inevitably implies a combination of 
direct and indirect studies, in which cosmological probes play an important role. Such probes provide 
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supplementary tools to speeify the expeeted properties of SUSY partieles both for their effieient seareh 
at aeeelerators and for their ability to explain eosmologieal dark matter. In the latter ease the existenee 
of new eonservation laws in SUSY models, sueh as R-parity, leads to the stability of the lightest partiele 
that possesses the eorresponding eonserved property. Created in the early Universe sueh partieles should 
survive to the present day and play the role of dark-matter eandidates. This is the simplest and most 
widely known eosmologieal impaet of Supersymmetry. 

However the fundamental basis of supersymmetrie models eontains a mueh wider field of 
eosmologieal eonsequenees. This basis ean naturally inelude physieal meehanisms for inflation and 
baryosynthesis, and eould beeome the physieal grounds for modem inflationary eosmology with 
baryosynthesis and dark matter. Sueh meehanisms imply physies of the very early Universe, whieh 
ean hardly be probed by direet experimental means, henee the neeessity to extend the set of theoretieal 
tools to study these phenomena. Moreover the absenee of positive results in the experimental seareh for 
SUSY partieles at aeeelerators ean move their mass seale to values so high that only indireet methods of 
study would be appropriate. 

Here, following the general framework of diseussion of eosmologieal effeets [1,2], we eonsider 
some nontrivial features of sueh effeets in their possible applieation to supersymmetrie models, paying 
speeial attention to predietions of nonequilibrium partieles and primordial nonlinear struetures, in whieh 
Primordial Blaek Holes (PBH), effeets of PBH evaporation and even antimatter domains in the baryon 
asymmetrieal Universe ean provide nontrivial probes for Supersymmetry in astronomieal observations 
and spaee experiments. 

After a brief review of the eosmologieal traees of partiele models and of their possible relationship 
with SUSY (Seetion 2), we turn to PBHs that ean play the role of a unique theoretieal tool in studies 
of physies of the very early Universe. We eonsider some meehanisms of PBH formation (Seetion 3) 
and possible effeets of SUSY partieles in PBH evaporation (Seetion 4). Metastable SUSY partieles 
or evaporating PBHs are the souree of nonequilibrium partieles after Big Bang Nueleosynthesis that 
ean influenee the primordial ehemieal eomposition or provide a non-thermal produetion of dark matter 
partieles (Seetion 5). A sueeession of phase transitions in the inflationary Universe ean result in 
primordial nonlinear struetures (Seetion 6). A profound signature of a nonhomogeneous SUSY-based 
meehanism for baryosynthesis is the appearanee of antimatter domains in a baryon asymmetrie Universe 
(Seetion 7). We eonelude that the development of supersymmetrie models perfeetly illustrates the might 
of the basie prineiples of eosmopartiele physios [3-8], whieh studies the physieal basis of the modem 
eosmology and oombines probes from eosmology, astrophysios and experimental physios in a systematio 
way (Seetion 8). 

2. Cosmological traces of particle models 

Let’s start following [1,2] (see also [9-13]) with a general diseussion of a variety of eosmologieal 
eonsequenees of partiele models and of their possible implementation in the SUSY ease. 

If R-parity is eonserved, the Lightest Supersymmetrie Partiele (LSP) should be stable. If these 
partieles were oreated at the early stages of eosmologieal evolution, they should survive in the Universe 
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until the present time. Therefore the simplest eosmologieal effeet of supersymmetrie models is a gas of 
new stable massive partieles, originated from the early Universe. 

For partieles with mass m, at high temperature ^ T > m the equilibrium eondition, n ■ av ■ t > 1 
is valid, if their annihilation eross seetion a > l/{mmpi) is suffieiently large to establish equilibrium. 
At T < m sueh partieles go out of equilibrium and their relative eoneentration freezes out. Weakly 
interaeting speeies deeouple from the plasma and the radiation at T > m, when n ■ av ■ t ^ 1, i.e. at 
Tdec ~ This is the main idea of the ealeulation of primordial abundanees for WIMP-like dark 

matter eandidates sueh as the neutralino (see e.g. [5,8] for details). 

The maximal temperature whieh is reaehed in an inflationary Universe is the reheating temperature, 
Tr, after inflation. So, very weakly interaeting partieles with annihilation eross seetion a < l/{Trmpi), 
as well as very heavy partieles with mass m ^ T^. eannot be in thermal equilibrium, and the detailed 
meehanism of their produetion should be eonsidered to ealeulate their primordial abundanees. This is 
the ease for the gravitino, whieh has a semigravitational interaetion eross seetion. 

By definition, primordial stable partieles survive to the present day and should be present in the 
modem Universe. The net effeet of their existenee is given by their eontribution to the total eosmologieal 
density. If they dominate at the matter dominated stage they play the role of eosmologieal dark matter 
that formed the large seale stmeture of the Universe. To be a dark matter eandidate, stable partieles 
should deeouple from the plasma and the radiation before the beginning of matter-dominated (MD) 
stage. This ean be satisfied for a very wide range of eross seetions of partiele interaetion with matter 
- from superweak to strong. The miraele of Weakly Interaeting Massive partieles (sueh as neutralino) 
is that their annihilation eross seetion provides the frozen out eoneentration that explains the observed 
dark matter density. For the eross seetions typieal of weak interaetion, a direet seareh for primordial 
dark-matter partieles in underground experiments is possible (see [14-17] and referenees therein). 

Dark-matter partieles ean annihilate in the Galaxy, eontributing by their annihilation produets to the 
flux of eosmie rays and gamma radiation. It was first found in [18] that sueh an effeet of dark matter 
annihilation ean provide a sensitive tool in indireet studies of dark matter by the measurements of eosmie 
positron and gamma baekground fluxes. 

There may be several types of new stable partieles. Then multi-eomponent dark matter seenarios are 
realized. It is interesting that sueh a multi-eomponent seenario ean naturally arise in the heterotie-string 
phenomenology, whieh ean embed both LSP and other possible types of stable partieles or objeets. 

The meehanism of eompaetifieation and symmetry breaking leads to the predietion of homotopieally 
stable objeets [19] and shadow matter [20], giving a wide variety of new types of dark matter eandidates. 

The embedding of the symmetry of the Standard model and Supersymmetry within the heterotie string 
phenomenology ean be also aeeompanied by the predietion of a fourth generation of quarks and leptons 
[21] with a stable massive 4th neutrino [22] and possibly a stable quark. 

Indeed, the eomparison between the rank of the unifying group Eq (r = 6) and the rank of the 
Standard model (r = 4) implies the existenee of new eonserved eharges and new (possibly striet) gauge 
symmetries. A new striet gauge U(l) symmetry (similar to the U(l) symmetry of eleetrodynamies) is 


1 


Throughout the paper we use the units ?i = c = /c = 1, if it is not specified otherwise . 
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possible, if it is ascribed to the fermions of the 4th generation. This hypothesis explains, in particular, 
the difference between the three known types of neutrinos and the neutrino of the 4th generation. The 
latter possesses a new gauge charge and, being a Dirac particle, cannot have a small Majorana mass due 
to see-saw mechanism. If the 4th neutrino is the lightest particle of the 4th quark-lepton family, the strict 
conservation of the new charge makes it absolutely stable. Following this hypothesis [21] the quarks and 
the leptons of the 4th generation are the source of a new long-range interaction (^/-electromagnetism), 
similar to the electromagnetic interaction of ordinary charged particles. 

New particles with an electric charge and/or a strong interaction can form anomalous atoms and 
be contained in the ordinary matter as anomalous isotopes. For example, if the lightest quark of 4th 
generation is stable, it can form stable charged hadrons, serving as nuclei of anomalous atoms of 
e.g. anomalous hydrogen and helium [23-28]. The experimental upper limits on the concentration 
of anomalous isotopes and especially of anomalous hydrogen severely constrains the possibility of new 
stable charged particles and practically rules out such particles with charges -i-l or -1. However, particles 
with charge -2 can exist and avoid these constraints and be hidden in 0-helium atoms - their neutral 
atom-like bound states with primordial helium nuclei. The 0-helium hypothesis can explain some 
puzzles of dark matter searches, challenging experimental searches for stable doubly charged particles 
at the LHC [1,29]. The -2 charge component of 0-helium can originate from stable anti-U quarks of a 
4th generation, which can form a A-like (UUU) state bound by chromoelectric forces [27]. 

Therefore together with SUSY candidates some other types of dark matter candidates - massive stable 
4th generation neutrinos, as well as nuclear-interacting 0-helium dark atoms, made of stable (anti-)U 
quarks of the 4th generation, are possible in the heterotic string phenomenology, the framework of which 
favors a multicomponent analysis of dark-matter effects. 

In the multi-component dark matter scenario a detailed study of the distribution of particles in space, 
of their condensation in galaxies, of their capture by stars. Sun and Earth, as well as of the effects of 
their interaction with matter and of their annihilation, provides a sensitive probe for the existence of even 
subdominant components. In particular, though stable neutrinos of the 4th generation with masses about 
50 GeV are predicted to be the subdominant form of dark matter, contributing less than 0,1 % to the 
total dark matter density [18,22,30-36], this possibility can be ruled out by direct experimental searches 
for WIMPs (see [14-17] and references therein) and by studies of the effects of the annihilation of 4th 
generation neutrinos and antineutrinos in the Galaxy in the galactic gamma-background. 4th generation 
neutrino annihilation inside the Earth should also lead to a flux of underground monochromatic neutrinos 
of known types, which can be excluded by the analysis of the existing data from underground neutrino 
detectors. 

Primordial unstable particles with a lifetime shorter than the age of the Universe, r < tjj, cannot 
survive to the present day. But, if their lifetime is sufficiently large to satisfy the condition r S> {iripi / m) ■ 
(1/m), their existence in the early Universe can lead to direct or indirect traces. The cosmological flux 
of decay products contributing to the cosmic- and gamma-ray backgrounds represents the direct trace of 
unstable particles. If the decay products do not survive to the present time their interaction with matter 
and radiation can cause indirect effects in the light element abundance or in the fluctuations of thermal 
radiation. 
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If the particle lifetime is much shorter than 1 s multi-step indirect traces are possible, provided that 
particles dominate in the Universe before their decay. During the dust-like stage of their dominance 
black hole formation takes place, and the spectrum of such primordial black holes traces the particle 
properties (mass, frozen concentration, lifetime) [37]. Particle decay in the end of the dust-like stage 
influences the baryon asymmetry of the Universe. In any case logical chains within a given cosmic 
phenomenology link the predicted properties of even unstable new particles to the effects accessible in 
astronomical observations. Such effects may be important for the analysis of the observational data. 

The parameters of new stable and metastable particles are also determined by the pattern of symmetry 
breaking. This pattern is reflected in a succession of phase transitions in the early Universe. First 
order phase transitions proceed through bubble nucleation, which can result in black hole formation 
[38] (see e.g. [39] and [40] for review and references). Phase transitions of the second order can lead 
to the formation of topological defects, such as walls, strings or monopoles. The observational data 
put severe constraints on magnetic monopoles [41,42] and cosmic wall production [43], as well as on 
the parameters of cosmic strings [44,45]. The structure of the cosmological defects can be changed in 
successive phase transitions. More complicated forms such as walls-surrounded-by-strings can appear. 
Such structures can be unstable, but their existence can leave a trace in the nonhomogeneous distribution 
of dark matter and give rise to large scale structures of nonhomogeneous dark matter such as archioles 
[46-48]. Primordial Black Holes, whose hypothetical existence was first formulated by Zeldovich and 
Novikov [49], represent a profound signature of such structures. 

3. The reflection of high-energy physics in the PBH spectrum 

Primordial Black Holes (PBHs) are a very sensitive cosmological probe of physics phenomena 
occurring in the early Universe. They could be formed by many different mechanisms (see e.g. [1,2] 
for a review and reference). 

After formation PBHs should remain in the Universe and, if they survive to the present day, they 
should represent a specific form of dark matter (see e.g. [50]). PBH evaporation by Hawking radiation 
[51] makes them a source of products of evaporation, which contain any type of particles that can exist in 
our space-time, both known and unknown. For a wide range of parameters the predicted effect of PBHs 
contradicts the data and puts restrictions on the mechanisms of PBH formation and on the underlying 
physics of the very early Universe. On the other hand, for some fixed values of their parameters, PBHs 
or their evaporation can provide a nontrivial solution to some astrophysical problems. 

Here we outline, following [1,2,50], the relationship of new physics with the mechanisms of PBH 
formation and the possible reflection of its parameters in the PBH spectrum. 

3.1. PBHs from superheavy metastable particles 

The formation of a black hole is highly improbable in a homogeneous expanding Universe, since it 
implies metric fluctuations of order 1. For Gaussian metric fluctuations with a dispersion (5^) 1, the 

probability of such fluctuations is determined by the exponentially small tail of the high-amplitude part 
of the distribution. This probability can even be more suppressed in the case of non-Gaussian fluctuations 
[ 2 ]. 
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If the Universe has an equation of state p = 7 e, with the numerieal factor 7 being in the range 
0 < 7 < 1 , the probability to form a black hole from fluctuations within the cosmological horizon is 
given by [52] 

lUpsH ocexp . ( 1 ) 

It provides an exponential sensitivity of the PBH spectrum to the softening of the equation of state in the 
early Universe (7 —)■ 0) or to the increase of the ultraviolet part of the spectrum of density fluctuations 
((5^) —>■ 1). These phenomena can appear as cosmological consequence of particle theory. 

It was first noticed in [53] that the dominance of superheavy metastable particles with lifetime r ^ 1 s 
in the Universe before their decay at f < r can result in the formation of PBHs, remaining in the 
Universe after the particles have decayed and keeping some information on the particle properties in 
their spectrum. 

After reheating, at T < Tq = rm, particles with mass m and relative abundance r = n/ur (where 
n is the frozen-out concentration of particles and Ur is the concentration of relativistic species) must 
dominate in the Universe before their decay. Dominance of these nonrelativistic particles at f > t^, 
where to = ^npi/T^, corresponds to a dust-like stage with an equation of state p = 0 , for which the 
particle density fluctuations grow as 

5it) = ^ oc f2/3 ( 2 ) 

P 

and the development of gravitational instability results in the formation of gravitationally bound systems, 
which decouple from the general cosmological expansion at 

t tf ^ (3) 

with S{tf) ~ 1 for fluctuations, which enter the horizon at t = ti>to with amplitude S{ti). 

The formation of these systems can result in black-hole formation either immediately after the 
system decouples from the expansion or as a result of the evolution of initially formed nonrelativistic 
gravitationally bound systems [54,55]. 

If the density fluctuation is especially homogeneous and isotropic, it directly collapses to PBH as 
soon as the amplitude of fluctuation grows to 1 and the system decouples from the expansion. The 
probability for direct PBH formation in the collapse of such homogeneous and isotropic configurations 
gives a minimal estimate of PBH formation at the dust-like stage. 

The mechanism [2,53] is effective for the formation of PBHs with masses in the interval 


Mo<M < M^hmax. (4) 

The minimal mass corresponds to the mass within the cosmological horizon at time t to, when 
particles start to dominate in the Universe and it is equal to [2,53] 

Mo = ^ptl ^ (5) 

3 rm 

The maximal mass is indirectly determined by the condition 


T = t{Mbhmax)S{Mbhmax) 


( 6 ) 
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that fluctuations at the considered scale M},hrnax, entering the horizon at t{Mi,hrnax) with an amplitude 
S{Mbhmax), can manage to grow to the nonlinear stage, decouple and collapse before particles decay at 
t = T. For a scale-invariant spectrum 5{M) = 5o the maximal mass is given by [2] 

Mbhmax = (7) 

tpi 

The direct mechanism of PBH formation can also be effective during a prereheating dust-like 
post-inflational stage of inflaton field oscillations [56]. 

3.2. Spikes from phase transitions during an inflationary stage 

A scale-dependent spectrum of fluctuations, in which the amplitude of small scale fluctuations is 
enhanced, can be another factor increasing the probability of PBH formation. The simplest functional 
form of such a spectrum is represented by a blue spectrum with a power law dispersion 

{5‘^{M)) oc (8) 

with an amplitude of fluctuation growing at small M for A; > 0. A realistic account of the existence of 
other scalar fields together with the inflaton during the period of inflation can give rise to spectra with 
distinctive scales, determined by the parameters of the considered fields and by their interaction. 

In a chaotic inflation scenario, the interaction of a Higgs field f with the inflaton rj can give rise to 
phase transitions during the inflationary stage, if this interaction induces a positive mass term +'^Tf(jfl. 
When in the course of slow rolling the amplitude of the inflaton decreases below a certain critical value 
Vc = the mass term in the Higgs potential 

Tn‘^ \ 7 /^ 

V{(l),v) = —( 9 ) 

changes sign and a phase transition takes place. Such a phase transitions during the inflationary stage 
lead to the appearance of characteristic spikes in the spectrum of initial density perturbations. These 
spike-like perturbations, at scales that cross the horizon from 60 to 1 e- folds before the end of inflation 
reenter the horizon during the radiation or dust-like era and could in principle collapse to form primordial 
black holes. The possibility of such spikes in a chaotic inflation scenario was first pointed out in [57] 
and implemented in [58] as a mechanism of PBH formation. 

If a phase transition takes place at e-folding N before the end of inflation, the spike re-enters the 
horizon at the radiation dominance (RD) stage and forms Black hole of mass 

9 

TTt 

M ^-f^exp{2Ar}, (10) 

-no 

where Hq is the Hubble constant during the period of inflation. 

If the spike re-enters the horizon during the matter dominance (MD) stage it should form black holes 
of mass 

~ ^exp{3Ar}. 

-no 


M 


( 11 ) 
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3.3. PBHs from a first-order phase transitions in the early Universe 

First-order phase transitions go through bubble nueleation. The simplest way to deseribe first order 
phase transitions with bubble ereation in the early Universe is based on a sealar field theory with two 
non degenerate vaeuum states. Being stable at the elassieal level, the false vaeuum state deeays due to 
quantum effeets, leading to a nueleation of bubbles of true vaeuum and their subsequent expansion [ 2 ]. 
The potential energy of the false vaeuum is eonverted into the kinetie energy of the bubble walls thus 
making them highly relativistie in a short time. The bubble expands till it eollides with another one. As 
it was shown in [39] a blaek hole may be ereated in a eollision of two bubbles. 

Just after the eollision, the mutual penetration of the walls up to a distanee eomparable with their 
width is aeeompanied by a signifieant inerease in potential energy [2]. Then the walls refleet eaeh other 
and aeeelerate baekwards. The spaee between them gets filled by the field in the false vaeuum state, 
eonverting the kinetie energy of the wall baek to the energy of the false vaeuum state and slowing down 
the veloeity of the walls. Meanwhile the outer area of the false vaeuum is absorbed by the outer wall, 
whieh expands and aeeelerates outwards. There is an instant when the eentral region of the false vaeuum 
is eausally separated fromm the walls [59] . If this false vaeuum bag shrinks in its oseillations within its 
gravitational radius, a blaek hole is formed. The mass of this PBH is given by (see [39]) 


Mbh — yiMbub 


( 12 ) 


where 71 ~ 10 “^ and M^ub is the mass that eould be eontained in the bubble volume at the epoeh of 
eollision for a full thermalization of bubbles. 

If inflation ends with a first-order phase transition, eollisions between bubbles of Hubble size in the 
pereolation regime lead to eopious PBH formation with masses 


Mo = = 


hor _ 7l ^pl 


2 H, 


end 


(13) 


where is the mass within the Hubble horizon at the end of inflation. Aeeording to ([39]) the initial 
mass fraetion of these PBHs is given by f3o ~ 71/e ~ 6 ■ 10“^. For example, for typieal value of 
Hend ~ 4 ■ the initial mass fraetion f3 eorresponds to PBHs with masses Mq ~ 1 g. 


4. Gravitino production by PBH evaporation 


Presently there is no observational evidenee proving the existenee of PBHs. However, even the 
absenee of PBHs provides a very sensitive theoretieal tool to study the physies of the early Universe. 
PBHs represent a nonrelativistie form of matter and their density deereases with the seale faetor a as 
a~^ oc T^, while the total density is cx: a~^ oc in the period of radiation dominanee (RD). As it must 
be formed within the horizon, a PBH of mass M ean be formed only after 


f(M) 


M _ M 

^pl — 9 “ • 

rupi 


(14) 
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If they are formed during the RD stage, the smaller are the masses of PBHs, the larger their relative 
eontribution to the total density during the modem MD stage. Therefore, even the modest eonstraint on 
the density of PBHs of mass M 


a,MM) = 


(15) 


pc 


in units of eritieal density pc = (SttG) from the eondition that their eontribution a{M) to the total 
density 


a(M) ^ 

ptot 

for ptot = Pc does not exceed that of dark matter 


^pbh{M) 


(16) 


cr(dT) — < tlptM — 0.23 


(17) 


eonverts into a severe eonstraint on their contribution 

^ Ppbh{M, tf) 
^ ~ Ptotitf) 


(18) 


during the period of their formation. If formed during the RD stage at = t{M), given by (14), which 
corresponds to the temperature Tf = rupi ^Jmpl/M, PBHs contribute to the total density at the end of the 
RD stage at tcq, corresponding to Tcq ~ leV, by a factor a{teq)/a{tf) = Tj/T^q = rUpi/Tgq \/rUpi/M 
larger than in the period of their formation. The constraint on f3{M), following from Eq.(17) is then 
given by 

/)(M) = a(M)^./^<0.23^./^. (19) 

mpi y TUpi TUpi y rUpi 

The possibility of PBH evaporation, revealed by S. Hawking [51], strongly influences the effects of 
PBHs. In the strong gravitational field near gravitational radius Vg of the PBH, the quantum effect of 
creation of particles with momentum p ~ l/r^ is possible. Due to this effect, the PBH becomes a 
black-body source of particles with temperature (in the units h = c = k = 1) 


T 


1 

SttGM 


^ lO^^GeV^. 

M 


The evaporation time scale can be written in the following form 


'TbH 


9*m^pi 


( 20 ) 


( 21 ) 


where is the number of effective massless degrees of freedom. For M < 10^^ g, it is less, than the age 
of the Universe. Such PBHs can not survive to the present day and the magnitude of Eq.(17) should be 
re-defined as the contribution to the total density at the moment of PBH evaporation. For PBHs formed 
during the RD stage and evaporated during the RD stage at f < tcq, the relationship Eq.(19) between 
/3(M) and a{M) is given by [37,60] 

( 22 ) 

The relationship between (3{M) and a{M) has a more complicated form if PBHs are formed during 
early dust-like stages [5,37,61,62], or if such stages take place after the PBH formation[5,62]. The 
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relative contribution of PBHs to the total density does not grow during the dust-like stage and the 
relationship between f3{M) and a{M) is model-dependent. A minimal model-independent factor 
a{M)/(3{M) follows from the enhancement taking place only during the RD stage between the first 
second of expansion and the end of RD stage at t^q, since radiation dominance in this period is supported 
by the light element abundances and by the spectrum of the Cosmic Microwave Background (CMB) 
[5,8,37,61,62] 

The effects of PBH evaporation make astrophysical data much more sensitive to the existence of 
PBHs. Constraining the abundance of primordial black holes can lead to invaluable information on 
cosmological processes, particularly as they are probably the only viable probe of the power spectrum at 
very small scales, which remain far from the CMB and from the Large Scale Structures (LSS) sensitivity 
ranges. To date, only PBHs with initial masses between ~ 10® g and ~ 10^® g have led to stringent 
limits (see e.g. [37,63-66]) from the consideration of the entropy per baryon, the deuterium destruction, 
the “^He destruction and the cosmic rays currently emitted by the Hawking process [51]. The existence 
of light PBHs should lead to important observable constraints, either through the direct effects of the 
evaporated particles (for initial masses between 10^^ g and 10^® g) or through the indirect effects of their 
interaction with matter and radiation in the early Universe (for PBH masses between 10® g and 10^"^ g). 

Several constraints on the density of PBHs have been derived in different mass ranges assuming 
the evaporation to standard model particles only: for 10® g < M < 10^^ g the entropy per baryon at 
nucleosynthesis was used [67] to obtain [3 < (10® g/M), for 10® g < M < 10^^ g the production of nn 
pairs at nucleosynthesis was used [68] to obtain (3 < 3 x 10“^^(10® g/M)^/®, for 10^® g < M < 10^^ g 
deuterium destruction was used [69] to obtain < 3 X 10 ^^(M/10^® for 10^^ g < M < 10^® g 
spallation of ^He was used [62,70] to obtain (3 < 3x 10“^^(M/10® g)®/^, for M 5 x 10^^ g the gamma 
and cosmic rays were used [71,72] to obtain f3 < 10“^®. Slightly more stringent limits were obtained in 
[73], leading io (3 < 10“^® for masses between 10® g and 10^® g and in [74], leading io j3 < 10“®® for 
M = 5 X 10^^ g. Gamma rays and antiprotons were also recently re-analyzed in [75] and [76], improving 
a little the previous estimates. Such constraints, related to phenomena occurring after nucleosynthesis, 
apply only for black holes with initial masses above ~ 10® g. Below this value, the only limits for a long 
time was the very weak entropy constraint (related with the photon-to-baryon ratio). 

However, since the evaporation products are created by the gravitational field, any quantum with 
a mass lower than the black hole temperature should be emitted, independently of the strength of its 
interaction. This could provide a copious production of superweakly interacting particles that cannot 
be in equilibrium with the hot plasma of the very early Universe. It makes evaporating PBHs a unique 
source of all the species which can exist in the Universe. 

To derive a limit in the initial mass range rupi < M < 10^^ g, gravitinos emitted by black holes® were 
considered in [78]. Gravitinos are expected to be present in all local supersymmetric models, which are 
regarded as the more natural extensions of the standard model of high energy physics (see, e.g., [79] for 
an introductory review). 


2 


The limits on the abundance of PBHs from gravitino production were first calculated in [77]. 



Symmetry 2015 , 7 


825 


Following [5,8,62,80] and [77,81] (but in a different framework and using more stringent eonstraints), 
limits on the mass fraetion of blaek holes at the time of their formation (/3 = ppbh/P tot) were derived in 
[78] using the produetion of gravitinos during the evaporation proeess. Depending on whether gravitinos 
are expeeted to be stable or metastable, limits are obtained using the requirement that they do not 
overelose the Universe or that the formation of light nuelei by the interaetions of ^He nuelei with a 
nonequilibrium flux of D,T,^He and ^He does not eontradiet the observations. This approaeh is more 
eonstraining than the usual study of photo-dissoeiation indueed by photons-photinos pairs emitted by 
deeaying gravitinos. It opened a new window for the upper limits on (5 below 10® g. 


5. Nonequilibrium particles 


Superweakly interacting gravitinos could not be in thermal equilibrium and thus the mechanisms of 
their production gives rise to appearance of nonequilibrium particles during the radiation dominated 
stage. Metastable particles, decaying via hadronic or electromagnetic channels after Big Bang 
nucleosynthesis, result in fluxes of photons and electron-positron pairs with energies much larger than 
the thermal energies in this period. Nucleon-antinucleon pairs from such decays represent a profound 
example of the creation of species that cannot be in equilibrium in the considered period. The interaction 
of primordial nuclei with nonequilibrium particles from electromagnetic and hadronic cascades leads 
to recoil nuclei and nuclear fragments with energies above the Coulomb barrier, and makes successive 
nuclear reactions possible. 

In the framework of minimal Supergravity (mSUGRA), the gravitino mass is, by construction, 
expected to lie around the electroweak scale, i.e. in the 100 GeV range. In this case, the gravitino is 
metastable and decays after nucleosynthesis, leading to important modifications of the nucleosynthesis 
paradigm. Instead of using the usual photon-photino decay channel, the study of [78] relied on the 
more sensitive gluon-gluino channel. Based on [82-86], the antiprotons produced by the fragmentation 
of gluons emitted by decaying gravitinos were considered as a source of nonequilibrium light nuclei 
resulting from collisions of those antiprotons on equilibrium nuclei. Then, ®Li, ^Li and ^Be nuclei 
production by the interactions of the nonequilibrium nuclear flux with “^He nuclei in equilibrium was 
taken into account and compared with data (this approach is supported by several recent analyses [87,88] 
which lead to similar results). The resulting Monte-Carlo estimates [85] lead to the following constraint 
on the concentration of gravitinos: 77 , 3/2 < 1.1 x 10“^^777^/^/^, where 7773/2 is the gravitino mass in GeV. 
This constraint has been successfully used to derive an upper limit on the reheating temperature [85]: 
Tp < 3.8 X 10® GeV. The consequences of this limit on cosmic-rays emitted by PBHs was considered, 
e.g., in [89]. In the approach of [78] this stringent constraint on the gravitino abundance was related 
to the density of PBHs through direct gravitino emission. The usual Hawking formula [51] was used 
for the number of particles of type i emitted per unit of time t and per unit of energy Q. Introducing 
the temperature defined by Eq. (20) T = hc^/{16'K‘^kGM) ^ (10^^g)/M GeV, taking the relativistic 
approximation for T^, and integrating over time and energy, the total number of quanta of type i can be 
estimated as: 


27 X 10^4 dT r x^dx 

647i^asuGRA Jt, Jm/T - (- 1 )^ 


( 23 ) 
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where T is in GeV, rupi ^ 10“^ g, x = Q/T, m is the partiele mass and asuGRA aeeounts for the 
number of degrees of freedom through M'^dM = —asuGRAdt where M is the blaek hole mass. Onee 
the PBH temperature is higher than the gravitino mass, gravitinos will be emitted with a weight related 
to their number of degrees of freedom. Computing the number of emitted gravitinos as a funetion of 
the PBH initial mass and matehing it with the limit on the gravitino density imposed by nonequilibrium 
nueleosynthesis of light elements leads to an upper limit on the PBH number density. If PBHs are 
formed during a radiation dominated stage, this limit ean easily be eonverted into an upper limit on fd 
by evaluating the energy density of the radiation at the formation epoeh. The resulting limit is shown on 
Fig. 1 and leads to an important improvement over previous limits, nearly independently of the gravitino 
mass in the interesting range. This opens a new window on the very small seales in the early Universe. 

It is also possible to eonsider limits arising in Gauge-Mediated-Susy-Breaking (GMSB) models [90]. 
Those alternative seenarios, ineorporating a natural suppression of the rate of flavor-ehanging neutral 
eurrents due to the low energy seale, prediet the gravitino to be the Lightest Supersymmetrie Partiele 
(LSP). The LSP is stable if R-parity is eonserved. In this ease, a limit was obtained [78] by requiring 
^ 3 / 2,0 < i-e. by requiring that the eurrent gravitino density does not exeeed the matter density. It 
ean easily be derived from the previous method, by taking into aeeount the dilution of gravitinos in the 
period of PBH evaporation and the eonservation of gravitino to keep speeifie entropy ratio, that [78]: 


/ 9 < 


n 


M,0 


JV3/2^VU 



1 

2 


(24) 


where N 3/2 is the total number of gravitinos emitted by a PBH with initial mass M, t^q is the end of the 
RD stage and tf = max{tform, tend) is teh time at whieh a non-trivial equation of state for the period 
of PBH formation is eonsidered, e.g. a dust-like phase whieh ends at tend [61]. The limit (24) does 
not imply the thermal equilibrium of the relativistie plasma in the period before PBH evaporation and is 
valid even for low reheating temperatures provided that the equation of state during the preheating stage 
is elose to relativistie. With the present matter density Um,o ~ 0.30 [91] this leads to the limit shown 
in Fig. 2 for 7713/2 = 10 GeV. Following (24) this limit seales with gravitino mass as Models of 
gravitino dark matter with 03 / 2,0 = ^gdm,o, corresponding to the case of equality in the above formula, 
were considered in [92,93]. 


6. Massive Primordial Black Holes from the collapse of closed walls 

A wide class of particle models possesses a symmetry breaking pattern, which can be effectively 
described by pseudo-Nambu-Goldstone (PNG) fields and which corresponds to the formation of an 
unstable topological defect structure in the early Universe (see [40] for a review and references). 
The Nambu-Goldstone nature of such an effective description reflects the spontaneous breaking of 
a global U(l) symmetry, resulting in a continuous degeneracy of vacua. The explicit symmetry 
breaking at a smaller energy scale changes this continuous degeneracy to a discrete vacuum degeneracy. 
The characteristics of the formed structures is different for phase transitions taking place during 
post-inflationary or inflationary stages. 
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Figure 1. Constraints of [78] on the fraction of the Universe going into PBHs (adapted from 
[37,63-65]). The two curves obtained with gravitino emission in mSUGRA correspond to 
^^ 3/2 =100 GeV (lower curve in the high mass range) and m 3/2 = 1 TeV (upper curve in the 
high mass range) 
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Figure 2. Constraints of [78] on the fraction of the Universe going into PBHs (adapted 
from [37,63-65]). The curve obtained with gravitinos emission in GMSB correspond to 
^^3/2 = 10 GeV and scales with gravitino mass as oc 
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6.1. Structures from a succession ofU(l) phase transitions 

At high temperatures such a symmetry breaking pattern implies a succession of second order phase 
transitions. In the first transition, continuous degeneracy of vacua leads, at scales exceeding the 
correlation length, to the formation of topological defects in the form of a string network; in the second 
phase transition, continuous transitions in space between degenerate vacua form surfaces: domain walls 
surrounded by strings. This last structure is unstable, but, as was shown in the example of the invisible 
axion [46-48], it is reflected in the large-scale inhomogeneity of the distribution of energy density of 
coherent PNG (axion) field oscillations. This energy density is proportional to the initial value of the 
phase, which acquires the dynamical meaning of the amplitude of the axion field, when the axion mass 
is switched on in the result of the second phase transition. 

The value of the phase changes by 27r around the string. This strong nonhomogeneity of phase leads 
to a corresponding nonhomogeneity of the energy density of coherent PNG (axion) field oscillations. 
The usual argument (see e.g. [94] and references therein) is essential only on scales corresponding to the 
mean distance between strings. This distance is small, being of the order of the scale of the cosmological 
horizon in the period during which PNG field oscillations start. However, since the nonhomogeneity of 
the phase follows the pattern of the axion string network, this argument misses large-scale correlations 
in the distribution of the energy density of oscillations’. 

Indeed, a numerical analysis of string networks (see review in [95]) indicates that large string loops 
are strongly suppressed and that a fraction of about 80% of string length, corresponding to long loops, 
remains virtually the same at all large scales. This property is the other side of the well known scale 
invariant character of string networks. Therefore the correlations in the energy density should persist on 
large scales, as was shown in [46-48]. 

The large-scale correlations in topological defects and their imprints in primordial inhomogeneities 
constitute an indirect effect of inflation, if phase transitions take place after the reheating of the Universe. 
Inflation provides in this case identical conditions for the phase transitions which takes place in causally 
disconnected regions. 

If phase transitions take place during the inflationary stage, new forms of primordial large scale 
correlations appear. The value of the phase after the first phase transition is inflated over the region 
corresponding to the period of inflation, while during inflation fluctuations of this phase change its 
initial value within regions of smaller size. Owing to such fluctuations, for a fixed value of 6 qq in the 
period of inflation with e-folding N = Qt) corresponding to the part of the Universe within the modem 
cosmological horizon, strong deviations from this value appear at smaller scales, corresponding to later 
periods of inflation with N <60. If < vr, the fluctuations can move the value of 6 * 7 v to 6 *Ar > vr in 
some regions of the Universe. After reheating as a result of the second phase transition these regions, 
which correspond to a vacuum with O^ac = 27r, are surrounded by the bulk of the volume with vacuum 
dvac = 0. As a result, massive walls are formed at the border between the two vacua. Since regions with 
dvac = 27r are confined, the domain walls are closed. As they fully enter the horizon, closed walls can 
collapse into black holes. 

This mechanism can lead to the formation of primordial black holes of arbitrarily large mass (up to 
the mass of AGNs [96], see for latest review [97]). Such black holes appear in the form of primordial 
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black hole clusters, exhibiting a fractal distribution in space [40,98,99]. Theyt can shed new light on the 
problem of galaxy formation [40,100]. 

6.2. Formation of closed walls in an inflationary Universe 

To describe a mechanism for the appearance of massive walls of a size essentially greater than the 
horizon at the end of inflation, let us consider a complex scalar field with the potential [40,96,98,99] 

V{^) = \M^-f/2f + 5V{e), (25) 

where = re*®. This field coexists with an inflaton field which drives the Hubble constant H during the 
inflational stage. The term 

51/(0) =A^(1-cos 0), (26) 

reflecting the contribution of instanton effects to the Lagrangian renormalization (see for example [101]), 
is negligible during the inflationary stage and during some period in the FLRW expansion. The omitted 
term (26) becomes significant when the temperature falls below the values T ~ A. The mass of the 
radial field component r is assumed to be sufficiently large with respect to H, which means that the 
complex field is in the ground state even before the end of inflation. Since the term (26) is negligible 
during inflation, the field has the form p) f /V2 ■ e*® and the quantity f9 acquires the meaning of a 
massless field. 

At the same time, the well established behavior of quantum field fluctuations on de Sitter background 
[102] implies that the wavelength of a vacuum fluctuation of every scalar field grows exponentially, 
having a fixed amplitude. Namely, when the wavelength of a particular fluctuation, in the inflating 
Universe, becomes greater than the average amplitude of this fluctuation freezes out at some 
non-zero value because of the large friction term in the equation of motion of the scalar field, whereas 
its wavelength grows exponentially. Such a frozen fluctuation is equivalent to the appearance of a 
classical field that does not vanish after averaging over macroscopic space intervals. Because the vacuum 
must contain fluctuations at every wavelength, inflation leads to the creation of more and more new 
regions containing a classical field of specific amplitude with a scale greater than H~^. In the case of 
an effectively massless Nambu-Goldstone field considered here, the averaged amplitude of the phase 
fluctuations generated during each e-fold (time interval H~fl is given by 

66 = H/27rf. (27) 

Let us assume that the part of the Universe observed inside the contemporary horizon = 

3000h“^Mpc was inflating, over Nu ~ 60 e-folds, from a single causally connected domain of size 
for which the average value of the phase is 9q. When inflation begins in this region, after one 
e-fold, the volume of the Universe increases by a factor . The typical wavelength of the fluctuation 69 
generated during every e-fold is equal to H~^. Thus, the whole domain H~^, containing 9o, after the first 
e-fold effectively becomes divided into separate, causally disconnected domains of size H~^. Each 
domain corresponds to an almost homogeneous phase value Oq ± 69. Thereby, more and more domains 
appear with time, in which the phase differs significantly from the initial value 9q. A crucial point is the 
appearance of domains with a phase 0 > vr. Appearing only after a certain period of time during which 
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the Universe exhibited exponential expansion, these domains turn out to be surrounded by a spaee with 
phase 6 < TT. The eoexistenee of domains with phases 6 < n and 6 > tt leads, in the following, to the 
formation of a large-seale strueture of topologieal defeets. 

The potential (25) possesses a (7(1) symmetry whieh is spontaneously broken at least after some 
period of inflation. Note that the phase fluetuations during the first e-folds may, generally speaking, 
transform into fluetuations of the eosmie mierowave radiation, and this will lead to restrictions on the 
scaling parameter /. This difficulty can be avoided by taking into account the interaction of the field 
99 with the inflaton field (i.e. by making parameter / a variable [40]). This spontaneous breakdown is 
holding through the condition of smallness of the radial mass, > H. At the same time the 

condition 

2/2 

me = — H (28) 

on the angular mass provides the freezing out of the phase distribution until some moment of the FRW 
epoch. After the violation of condition (28) the term (26) contributes significantly to the potential (25) 
and explicitly breaks the continuous symmetry along the angular direction. Thus, the potential (25) 
eventually has a number of discrete degenerate minima in the angular direction at the points 6 'mm = 
0, ±27r, idvr, ... . 

As soon as the angular mass m^ is of the order of the Hubble rate, the phase starts oscillating about 
the potential minimum, with different initial values being in various space domains. Moreover, in the 
domains with the initial phase tt < 6 < 2tt, the oscillations proceed around the potential minimum at 
Qmin = 27r, whereas the phase in the surrounding space tends to a minimum at the point 9min = 0. 
At the end of the decaying phase oscillations, the system contains domains characterized by the phase 
9min = 27r surrounded by space with 9min = 0. Apparently, if we move in any direction from the inside 
to the outside of the domain, we will unavoidably pass through a point where 9 = tt because the phase 
varies continuously. This implies that a closed surface characterized by the phase 9^aii = tt must exist. 
The size of this surface depends on the moment of the domain formation in the inflation period, while the 
shape of the surface may be arbitrary. The key point for the subsequent considerations is that the surface 
is closed. After reheating of the Universe, the evolution of domains with the phase 9 > tt proceeds on the 
background of the Friedman expansion and is described by the relativistic equation of state. When the 
temperature falls down to T* ~ A, an equilibrium state between the "vacuum" phase 9yac = 27r inside the 
domain and the 9yac = 0 phase outside of it is established. Since the equation of motion corresponding to 
potential (26) admits a kink-like solution (see [95] and references therein), which interpolates between 
two adjacent vacua 9yac = 0 and 9yac = 27r, a closed wall corresponding to the transition region at 0 = tt 
is formed. The surface energy density of a wall of width ~ 1/m ~ //A^ is of the order ^ of /A". 

Note that if the coherent phase oscillations do not decay for a long time, their energy density can play 
the role of CDM. This is the case, for example, in the cosmology of the invisible axion (see [94] and 
references therein). 

It is clear that immediately after the end of inflation, the size of domains which contain a phase 
9vac > 2.TT exceeds by far the horizon size. This situation is replicated in the size distribution of vacuum 


^ The existence of such domain walls in the theory of the invisible axion was first pointed out in [103]. 
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walls, which appear at the temperature ^ A whenee the angular mass mg starts to build up. Those 
walls, whieh are larger than the eosmologieal horizon, still follow the general FLRW expansion until the 
moment when they get eausally eonneeted as a whole; this happens as soon as the size of a wall beeomes 
equal to the horizon size Rh- Evidently, internal stresses developed in the wall after erossing the horizon 
initiate proeesses tending to minimize the wall surfaee. This implies that the wall tends, first, to aequire 
a spherieal shape and, seeond, to eontraet toward its eentre. For simplieity, we will eonsider below the 
motion of elosed spherieal walls 

The wall energy is proportional to its area at the instant of erossing the horizon. At the moment of 
maximum eontraetion, this energy is almost eompletely eonverted into kinetie energy [106]. Should the 
wall at the same moment be loealized within the gravitational radius, a PBH is formed. 

A detailed study of blaek hole (BH) formation was made in [96] . The results of these ealeulations are 
sensitive to the parameter A and to the initial phase 9u. As the A value deereases to ~ 1 GeV, still greater 
PBHs appear with masses of up to ~ 10^° g. A ehange in the initial phase leads to sharp variations in 
the total number of blaek holes. As was shown above, eaeh domain generates a family of subdomains in 
its elose vieinity. The total mass of sueh a eluster is only 1.5-2 times that of the largest initial blaek hole 
in this spaee region. Thus, the ealeulations eonfirm the possibility of formation of elusters of massive 
PBHs ( ~ IOOMq and above) in the earliest stages of the evolution of the Universe at a temperature of 
~ 1 — 10 GeV. These elusters represent stable energy density fluetuations around whieh baryonie matter 
(and eold dark matter) may eoneentrate in the subsequent stages, followed by the evolution into galaxies. 

It should be noted that additional energy density is supplied by elosed walls of small sizes. Indeed, 
beeause of the smallness of their gravitational radius, they do not eollapse into BHs. After several 
oseillations sueh walls disappear, leaving eoherent fluetuations of the PNG field. These fluetuations 
eontribute to a loeal energy density exeess, thus faeilitating the formation of galaxies. 

The mass range of formed PBHs is eonstrained by the fundamental parameters of the model / and 
A. The maximal BH mass is determined by the eondition that the wall does not dominate loeally before 
it enters the eosmologieal horizon. Otherwise, loeal wall dominanee leads to a superluminal a oc R 
expansion for the eorresponding region, separating it from the other parts of the Universe. This eondition 
eorresponds to the mass [40] 

(29) 

The minimal mass follows from the eondition that the gravitational radius of BH exeeeds the width of 
the wall and is equal to [40,98] 

Mrmn = f ■ (30) 

Closed wall eollapse leads to a primordial gravitational wave speetrum, peaked at [2] 

i/o = 3-10''(A//)Hz (31) 

with energy density up to 

Ugw ~ (32) 


4 


The motion of closed vacuum walls has been derived analytically in [104,105]. 
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At / lO^^GeV this primordial gravitational wave background can reach flew ~ 10 For the 
physically reasonable values of 

1 < A < 10®GeV (33) 

the maximum of the spectrum is at 

3 ■ 10-3 < I/O < 3 ■ lO^Hz. (34) 

Another profound signature of the considered scenario consists in gravitational wave signals from the 
merging of PBHs in the PBH cluster. These effects can provide tests of the considered approach in the 
cLISA experiment. 

7. Antimatter in a baryon-asymmetric Universe? 

Primordial strong inhomogeneities can also appear in the baryon charge distribution. The 
appearance of antibaryon domains in the baryon asymmetrical Universe, reflecting the inhomogeneity 
of baryosynthesis, is a profound signature of such a strong inhomogeneity [107]. In the model of 
spontaneous baryosynthesis (see Ref. [108] for a review) the possibility of the existence of antimatter 
domains, surviving to the present day in an inflationary Universe with inhomogeneous baryosynthesis 
was considered in [109]. 

The mechanism of spontaneous baryogenesis [110] implies the existence of a complex scalar field 
X = (//v^) exp (6*) carrying the baryonic charge. The U{1) symmetry, which corresponds to the 
baryon charge, is spontaneously and explicitly broken. The explicit breakdown of the U (1) symmetry is 
caused by the phase-dependent term 


1/(0) = A^(1-cos0). (35) 

The possible baryon- and lepton-number violating interaction of the field x with matter fields can have 
the following structure [108] 

^ = gxQL+ h..c., (36) 

where the fields Q and L represent a heavy quark and a lepton, coupled to the ordinary matter fields. 
In the early Universe, at a time when the friction term, induced by the Hubble constant, becomes 
comparable to the angular mass thq = the phase 9 starts to oscillate around the minima of the 
potential and decays into matter fields according to (36). The coupling (36) gives rise to the following 
[108]: as the phase starts to roll down in the clockwise direction (Fig. 3), it preferentially creates an 
excess of baryons over antibaryons, while the opposite is true as it starts to roll down in the opposite 
direction. 

The fate of such antimatter regions depends on their size. If their physical size is larger than the 
critical size Lc = Sh? kpc [109], they survive annihilation with surrounding matter. The evolution of 
sufficiently dense antimatter domains can lead to the formation of antimatter globular clusters [112]. 
The existence of such a cluster in the halo of our Galaxy should lead to the pollution of the galactic 
halo by antiprotons. Their annihilation can reproduce [113] the observed galactic gamma background in 
the range of tens to hundreds of MeV. The prediction of an antihelium component in cosmic rays [114], 
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Figure 3. The inflationary evolution of the phase (taken from [11 1]). The phase O^q sits in the 
range [tt, 0] at the beginning of inflation and makes Brownian step 56eff = Hinfi/{2rfeff) 
at eaeh e-fold. The typieal wavelength of the fluetuation 50 is equal to The whole 

domain containing the phase On gets divided, after one e-fold, into causally 

disconnected domains of radius Each new domain contains an almost homogeneous 

phase value 0n-i = On Or SOefj. Every successive e-fold, this process repeats in every 
domain. 

accessible to searches for cosmic ray antinuclei in the AMS02 experiment, as well as of antimatter 
meteorites [115] provides a direct experimental test for this hypothesis. The possibility of formation of 
dense antistars within an extension of the Affleck-Dine scenario of baryogenesis and the strategies for 
their search were recently considered in [1 16]. 

So primordial strong inhomogeneities in the distribution of total, dark matter and baryon density in 
the Universe are a new important phenomenon in cosmological models based on physics beyond the 
Standard model. 

8. Supersymmetry in the context of Cosmoparticle physics 

Observational cosmology offers strong evidence favoring the existence of new physics, and suggests 
experimental approaches to their investigation. Cosmoparticle physics [3-5,8], studying the physical, 
astrophysical and cosmological impact of new laws of Nature, explores new forms of matter and their 
physical properties. The development of SUSY models follows the main principles of Cosmoparticle 
physics and offers a great challenge for theoretical and experimental research. Physics of dark matter in 
all its aspects plays important role in this process. 

The necessity to extend the Standard model by supersymmetry has serious theoretical reasons. 
Aesthetically, because it helps to achieve full unification for the Standard model; practically, because it 
removes its internal tensions. Supersymmetry can also provide a complete physical basis for cosmology. 
It can pretend to explain inflation, baryosynthesis and nonbaryonic dark matter. The white spots in the 
representations of supersymmetric models correspond to new unknown particles. The extension of the 
symmetry of the gauge group introduces new gauge fields, mediating new interactions. Global symmetry 
breaking results in the existence of Goldstone boson fields. 

In particle physics, direct experimental probes for the predictions of particle theory are the most 
attractive. The predicted supersymmetric partners of known particles are accessible to experimental 
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search at accelerators if their masses are within a few-TeV range. However, the predictions related to 
higher energy scale need non-accelerator or indirect means for their test. 

Cosmoparticle physics offers complementary tools via indirect and non-accelerator direct searches for 
new physics and its possible properties. In experimental cosmoarcheology, data can be obtained, which 
links the new physics with astrophysical observations (see [117]). In experimental cosmoparticle physics 
the parameters, fixed from the consitency of cosmological models and observations, define the level at 
which the new types of particle processes should be searched for (see [118]). These basic principles of 
cosmoparticle physics have been widely implemented in the development of supersymmetric models. 

The theories of everything should provide a complete physical basis for cosmology and naturally 
involve supersymmetry. The problem is that string theory [119] is now in the form of "theoretical theory", 
for which many doubt that experimental probes exist, cosmoparticle physics can remove these doubts. 

For a long time, scenarios with Primordial Black holes belonged dominantly to cosmological 
"fantasies", as they provided restrictions on physics of the very early Universe from the contradiction 
of their predictions with observational data. Even this negative type of information makes PBHs an 
important theoretical tool. Being formed in the very early Universe as an initially nonrelativistic form 
of matter, PBHs should have increased their contribution to the total density during the RD stage 
of expansion, while the effect of PBH evaporation should have strongly increased the sensitivity of 
astrophysical data to their presence. Indeed, astrophysical constraints on hypothetical sources of cosmic 
or gamma rays, on light element abundances and on the spectrum ofthe CMB can be linked to restrictions 
on superheavy particles in the early Universe and on first- and second-order phase transitions, thus 
making astrophysical data a sensitive probe of particle symmetry structure and of the pattern of its 
breaking at superhigh energy scales. 

The gravitational mechanism of particle creation in PBH evaporation makes evaporating PBHs a 
unique source of any species of particles which can exist in our space-time. At least theoretically, PBHs 
can be treated as the source of such particles, which are strongly suppressed in any other astrophysical 
mechanism of particle production, either due to a very large mass of these species, or owing to their 
superweak interaction with ordinary matter. 

By construction, astrophysical constraints exclude effects predicted to be larger than observed. At the 
edge such constraints convert into an alternative mechanism for the observed phenomenon. At some fixed 
values of parameters, the PBH spectrum can play a positive role and shed new light on old astrophysical 
problems. 

Common sense dictates that PBHs should have small sub-stellar mass. The formation of PBHs within 
the cosmological horizon, which was very small in the very early Universe, seems to argue for this 
viewpoint. However, phase transitions during the inflationary stage can provide spikes in the spectrum 
of fluctuations at any scale, or lead to the formation of closed massive domain walls of any size. 

In the latter case, the existence of primordial clouds of massive black holes around an 
intermediate-mass or supermassive black hole is possible. Such clouds have a fractal spatial distribution. 
This approach suggests a radically new scenario for galaxy formation. Traditionally, the Big Bang 
model assumes a homogeneous distribution of matter at all scales, whereas the appearance of observed 
inhomogeneities is related to the growth of small initial density perturbations. However, the analysis 
of the cosmological consequences of the particle theory indicates the possible existence of strongly 
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inhomogeneous primordial structures in the distribution of both dark matter and baryons. These 
primordial structures represent a new factor in galaxy formation theory. Topological defects such 
as the cosmological walls and filaments, primordial black holes, archioles in the models of axionic 
CDM, and inhomogeneous baryosynthesis (leading to the formation of antimatter domains in a 
baryon-asymmetric Universe [5,8,40,107,109,110,112-115,120-125]) are an incomplete list of possible 
primary inhomogeneities inferred from the existing elementary particle models. 

We can conclude that, within the modern cosmological paradigm, from the very beginning to the 
present time, the evolution of the Universe was governed by physical laws which we still don’t fully 
know. These laws must come from a fundamental particle symmetry beyond the Standard model and they 
imply the use of methods of cosmoparticle physics for their study. Cosmoparticle physics originates from 
the well-established relation between the microscopic and the macroscopic descriptions in theoretical 
physics. This is reminiscent of the links between statistical physics and thermodynamics, or between 
electrodynamics and the theory of the electron. At the end of the XXth Century, a new instance of this 
kind of relationship was realized. It came both from the cosmological necessity to go beyond the world of 
known elementary particles to settle the physical grounds for inflationary cosmology with baryosynthesis 
and dark matter, and from the necessity for particle theory to use cosmological tests as an important and 
in many cases unique way to probe its predictions. The development of supersymmetric models perfectly 
reflects this direction of fundamental knowledge. 
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